Electric control of a {Fe^} single-molecule magnet in a single-electron transistor 
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Using first-principles methods we study theoretically the properties of an individual {Fe^} single- 
molecule magnet (SMM) attached to metallic leads in a single-electron transistor geometry. We show 
that the conductive leads do not affect the spin ordering and magnetic anisotropy of the neutral 
SMM. On the other hand, the leads have a strong effect on the anisotropy of the charged states of 
the molecule, which are probed in Coulomb blockade transport. Furthermore, we demonstrate that 
an external electric potential, modeling a gate electrode, can be used to manipulate the magnetic 
properties of the system. For a charged molecule, by localizing the extra charge with the gate 
voltage closer to the magnetic core, the anisotropy magnitude and spin ordering converges to the 
values found for the isolated {Fe^} SMM. We compare these findings with the results of recent 
quantum transport experiments in three-terminal devices. 



I. INTRODUCTION 

In recent years molecular spintronics has emerged 
as one of the most active areas of research within 
magnetism at the atomic scaled— Progress in the 
field is driven in part by advances in chemical de- 
sign and synthesis, which allow the realization of in- 
teresting magnetic molecules with desired electronic 
and magnetic properties. A second essential feature 
of ongoing research is the improved ability of in- 
tegrating individual magnetic molecules into solid- 
state nano-electronic devices. 

Typically magnetic molecules have long spin- 
relaxation times, which can be utilized in high- 
density information storage. They are also usually 
characterized by a weak hypcrfine interaction with 
the environment, resulting in long spin coherence 
times, which is an essential property for applications 
in quantum information processing. Single-molecule 
magnets (SMMs) are a special class of spin-ordered 
and/or magnetically active molecules characterized 
by a relatively high molecular spin and large mag- 
netic anisotropy energy^ The latter lifts the spin 
degeneracy even at zero magnetic field, and favors 
one particular alignment of the spin, making the 
molecule a nanoscale magnet. 

One of the goals of molecular spintronics is to 
address the magnetic states of individual magnetic 
molecules with electric fields and electric currents. 
In the last six years experimental efforts toward this 
goal have considered different classes of magnetic 
molecules and strategies to incorporate them into 
electric nano-circuits. A particularly interesting di- 
rection focuses on quantum transport in a single- 



electron transistor (SET), a three-terminal device 
where a SMM bridges the nanogap between two con- 
ducting nano-lcads, and can further be electrically 
manipulated by the gate voltage of a third nearby 
electrode. In the regime of weak coupling to the 
leads, the electric charge on the central SMM is 
quantized and can be controlled by the external gate. 
When Coulomb blockade is lifted by either gate or 
bias voltages, transport occurs via tunneling of sin- 
gle electrons in and out of the SMM. Therefore a 
study of transport in this geometry can provide de- 
tailed information on the magnetic properties of in- 
dividual SMMs, both when the molecule is neutral 
and when it is electrically charged. 

Early SET experiments on SMMs&2 focused on 
the archetypal SMM {Mni2j-acetat o 7 ' 10 , character- 
ized by the ground-state spin S = 10 and a large 
magnetic anisotropy barrier of approximately 50 K. 
Unfortunately these experiments and studies of self- 
assembled molecules on gold surfaces^ 1 - have shown 
that the magnetic properties of {Mn^} complexes 
are extremely fragile and easily disrupted when the 
molecule is attached to metallic leads or surfaces. 

More recently, another class of SMMs, namely the 
tctranuclear {Fe^} molecule, has emerged as a can- 
didate in molecular spintronics that does not suf- 
fer the drawbacks of {Mn^}. The properties of 
{Fe^} in its neutral state are well studied in the 
crystal phas o 12 ' 13 , and include a molecular ground- 
state spin S = 5 and an intermediate magnetic 
anisotropy barrier s=s 15 K. In contrast to what hap- 
pens with {M7112}, these magnetic characteristics 
remain stable when the molecule is deposited on a 
gold surface i 14 ' 15 Furthermore, its tripodal ligands 
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are shown to be advantageous for the preparation 
of single-molecule electronic devices. Indeed, recent 
three-terminal quantum transport experiments^— , 
with {Fe^} as the central island of a SET, show that 
this molecule behaves indeed as a nanoscale magnet, 
even when it is wired to metallic leads. The mag- 
netic anisotropy is significantly affected by adding 
or subtracting single charges to the molecule^, an 
operation that can be performed reversibly with the 
gate voltage. More refined techniques^ allow the 
extraction of the magnetic anisotropy of the neutral 
and charged molecule from the transport measure- 
ments with unprecedented accuracy. 

In this paper we carry out density functional the- 
ory (DFT) calculations to evaluate the magnetic 
properties of a {Fe^} connected to gold electrodes 
and under the effect of an external electric field rep- 
resenting a gate voltage. The geometry considered 
here is supposed to model the typical situation real- 
ized in current SET experiments, although some de- 
tails might be different. The main aim of our work 
is to investigate theoretically how the spin order- 
ing and the magnetic anisotropy of {Fe^} are af- 
fected by weak coupling to the leads, both when the 
molecule is in its neutral state and when a single 
charge is added to or subtracted from the device. A 
second important objective of the paper is a theo- 
retical analysis of how these magnetic properties can 
be modified and controlled by means of an external 
electric potential representing a gate electrode. 

Although a full-fledged first-principles study of 
quantum transport is beyond the scope of the 
present paper, as we explain below, we believe that 
our analysis of the charged states under the effect of 
an external electric field is useful to develop meth- 
ods to compute the tunneling conductance within 
a master equation formalism. Ref. [l9| introduced a 
DFT description of the neutral and charged states 
of an isolated {Mni 2 } SMM, which were then used 
in a master equation formalism for quantum trans- 
port. Here the coupling to the leads was treated with 
phenomcnological tunneling amplitudes taken form 
experiment. In Ref. I20I charge transport was studied 
by means of a non-equilibrium Green's function ap- 
proach, which included the presence of metallic leads 
in the regime of strong coupling. Transport calcula- 
tions carried out in the same transport regime but 
based on microscopic tight-binding models have also 
been considered^ This approach has the drawback 
that charging effects, essential for the description 
of SET experiments, cannot be adequately incorpo- 
rated. 

In this sense the present paper is a further contri- 
bution to these early attempts to use first-principles 
methods based on DFT to investigate quantum 



transport in a SET with a SMM. We are aware 
that the use of DFT can be problematic when it 
comes to describing the electronic structure and 
Fermi level alignment of molecules coupled to exter- 
nal electrodes, particularly when charged states are 
involved. Assessing the limitations of DFT in this 
context is also a technical aim of the present pa- 
per. In this regard, we discuss potential uncertain- 
ties in the relative level alignment of the electrode 
and molecular states. 

The main findings of our analysis are the follow- 
ing. The SMM {Fe^} in its neutral state is indeed 
quite robust against the presence of metallic leads: 
both spin ordering and magnetic anisotropy are es- 
sentially identical to the one of the isolated molecule. 
For the case of a charged molecule, the effect of the 
leads is more complex. Our calculations show that 
the extra charge tend to reside primarily on the lig- 
ands between molecule and metallic leads, and only 
minimally on the magnetic core. As a result the ad- 
dition of electrons affects the magnetism of molecule 
(specifically the magnetic anisotropy) considerably 
less than when the molecule is isolated. We find 
that an external gate voltage can be used to localize 
the extra charge closer to the central magnetic core, 
and in this case the magnetic characteristics of the 
device converge to the ones of the isolated {FtA] 
SMM. 

The organization of this paper is as follows. In 
Sec. [IT] we present an overview of the theoretical 
and computational approach employed in this work. 
The electronic and magnetic properties of different 
charge states of the isolated {Fe±} SMM is discussed 
in Sec. HID The effects of the metallic leads on the 
properties of the molecule are discussed in Sec. IIVI 
In Sec. [V] we discuss how a confining electric poten- 
tial affects the magnetic properties of the molecule. 
In Sec. lVII we compare the results of our calculations 
with recent SET transport experiments. Finally in 
Sec. IVIll we summarize our results. 



II. THEORETICAL BACKGROUND 

A. Spin Hamiltonian and the giant-spin model 

In a first approximation, the exchange interaction 
between the magnetic ions of a molecular magnet 
can be described by an isotropic Hciscnbcrg model 

H = ^ J iJ s i ■ s j > (!) 

where is the spin of the magnetic ion i and the 
constants Jjj describe the super-exchange coupling 



between ions i and j. Clearly, the validity of Eq. ([1} 
relies crucially on the assumption that each mag- 
netic ion is characterized by a well-defined quantum 
spin, localized at the ion position. There might be 
cases where the spin-polarization of the molecule is 
delocalized, where this assumption may break down. 
Such cases are especially probable when an excess 
tunneling electron is present. Once the exchange 
constants are known, the Hamiltonian can be diag- 
onalizcd. Since H is a sum of scalars in spin space, 
it commutes with the total spin S. Therefore the 
eigenvalues of S 2 and S z can be used to label the 
cigenstates of H. 

In the case of single-molecule magnets, the ground 
state (GS) of Eq. ([1]) is characterized by a relatively 
large spin S, and it is separated by a fairly large 
energy A J = Max{, 3 } [Jij] from excited states with 
different total spins. Thus, at low energies < A J the 
magnetic molecule behaves effectively as an atom 
with a relatively large spin 5, known as 'giant spin'. 
The approximation of restricting to the lowest spin 
multiplet is known as giant-spin model. According 
to Eq. ([1}, each spin multiplet is degenerate. In the 
next section we will discuss how spin-orbit interac- 
tion lifts this degeneracy, splitting the 25 + 1 states 
of the GS multiplet. 

Note that the ground-state spin S is not always 
the maximum value S'max allowed by Eq. (fTJ). Due 
to the presence of antifcrromagnctic components, 
the most common situation encountered in SMMs 
is an intermediate value 1 < S < S ma _ x , which quasi- 
classically corresponds to a ferrimagnetic spin con- 
figuration. Below we will show how all the physical 
quantities entering in the spin Hamiltonian of Eq. (fTJ 
and the value of the giant spin S can be calculated 
within DFT. 



B. Spin-orbit interaction and magnetic 
anisotropy barrier 

Spin-orbit interaction introduces terms to Eq. (fTJ 
that break rotational invariancc in spin space. Up 
to second-order perturbation theory, these terms, 
besides anisotropic corrections to the Heisenberg 
model, include the antisymmetric Dzyaloshinskii- 
Moriya spin exchange, and the single-ion magnetic 
anisotropy H ia = — ^(di • s i) 2 - Because of these 
terms, the total spin is no longer a good quantum 
number. Within the giant-spin model of SMMs, 
where the isotropic exchange is the dominant mag- 
netic interaction, the main effect of the spin-orbit 
interaction is to lift the spin degeneracy of the GS 
multiplet. To second-order perturbation theory, this 
can be described by the following anisotropy Hamil- 
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tonian for the giant spin operator S = (S x , S y , S z ) 

H = DSl + E{Sl - S 2 y ) . (2) 

The parameters D and E specify the axial and trans- 
verse magnetic anisotropy, respectively. If D < 
and \D\ >> \E\, which are defining properties for 
SMMs, the system exhibits an easy axis in the z- 
direction. In the absence of magnetic field, and ne- 
glecting the small transverse anisotropy term, the 
GS of Eq. ^ is doubly degenerate and it corre- 
sponds to the eigenstates of S z with eigenvalues ±5. 
To go from the state S z = +S to the state S z = — S 
the system has to surmount a magnetic anisotropy 
energy barrier AE = \D\S 2 . In addition, transitions 
which change the axial quantum numbers require 
some type of carrier to balance the change in spin. 
When the transverse term is not negligible and S z is 
not a good quantum number, we can still define an 
anisotropy barrier separating the two (degenerate) 
lowest energy levels as the energy difference between 
GS energy and the energy of the highest excited 
state. If D > the systems has a quasi-easy plane 
perpendicular to the z-axis without energy barrier. 

The anisotropy parameters D and E can be 
calculated within a self-consistent-field (SCF) one- 
particle theory (e.g. DFT or Hartree Fock), by in- 
cluding the contribution of the spin-orbit interac- 
tion. Here we summarize the main steps of the pro- 
cedure originally introduced in Ref. |22L (For more 
recent reviews see Refs. [23] andl24l.) 

The starting point are the matrix elements of the 
spin-orbit interaction (SOI) operator 

£/(r,p )S ) = -^s-px V$(r) (3) 

(p is the momentum operator; s is the electron 
spin operator; $ is the Coulomb potential and c 
is the speed of light), taken with respect to the 
unperturbed single-particle spinor wave functions 
IV'fco-) = \4>ka)\Xcr) , which are solutions of the SCF- 
approximation Schrodingcr equation 

H\^ ka ) =e ka \^ ka ) (4) 

Here ^>fc CT (r) = (r\\(f>k a ) is the orbital part of the 
wavefunction; the two-component spinors \xa), & = 
(1,2) are the eigenstates of sn, where the unit vector 
n = n(#, if) is an arbitrary quantization axis. 
The matrix elements can be written as^ 2 - 

Uko,k'a> = (tpka\U{r,p,s) \lp k >a>) 

= -i^ifikvlVt |<?W) (Xa\ s l \x*>) ,(5) 
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where the matrix elements of the operator V 
{V x ,V y , V z ) are defined as 



1 

2? 



v x \<t> Va .) 



! >ka 



dz 



d<t>, 



h'a 1 



dy 



Oka 



dy 



dz 



6) 



and cyclical. Note that Eq. ([6]) avoids the 
time-consuming calculation of the gradient of the 
Coulomb potential, replacing it with the calculation 
of the gradient of the basis functions in which <j>(v) 
is expanded. The above representation of the spin- 
orbit interaction arises by an integration by parts of 
the matrix element defined in Eq. (|6]). It is similar 
to the form of spin-orbit interaction that comes out 
of the Dirac equation. 

In the absence of an external magnetic field, the 
first-order perturbation-theory correction to the to- 
tal GS energy cause by the SOI is zero because of 
time-reversal symmetry. The second-order correc- 
tion is nonzero and can be written as 



be re-interpreted as the expectation values of the 
components of the giant-spin operator S for the spin- 
coherent state 1 5, n) with S = AN/2. 

Using the resolution of the identity in spin space, 
J2a Wo) (Xa\ = 1, we can write 



(Xi|si IX2) (X2\sj \xi) = 
(Xi I SiSj \xi) - (xi| Si |xi) (xi I Sj \Xi) 
(Si) (Sj) 



(xi\siSj \xi) 



(AiV) 2 



(11) 



and a similar expression for (x 2 | Sj Xi) (xi s j |X2 ) ■ 

With the help of Eqs. $TU\), (|TT|) and 
(Xo-I (si) 2 \x<r) = I/ 4 , Eq. © takes the form 



A 2 =a + J2'TiASi)(S j ) , 



(12) 



where a = J2ij(^u 2 + ^ii) 1S a constant indepen- 
dent of the quantization axis. The anisotropy tensor 
jij is given by 



A, 



where 



and 



EE 

i,j 



s r = (x*\si\x„>) 



(7) 



(8) 



E E 



^ka\Vi\(t)k'a')(4>k'a' \Vj \4>k<j) 



(9) 

where the sums over k and k' involve occupied and 
unoccupied states, respectively. 

Eq. (J7J is the central equation in the study of 
the magnetic anisotropy. Since the spin matrix el- 
ements sf a depend on the orientation of the arbi- 
trary axis of quantization n, so does also the total 
energy shift. This is precisely the origin of the mag- 
netic anisotropy brought about by SOI. 

We now consider the case of a closed-shell 
molecule, a system with a well defined HOMO- 
LUMO gap in order to avoid problems with partial 
occupancy, with AiV excess of majority spin elec- 
trons. 

Wc have the important relation 



(Xi\si \xi) = - (X2I Sj 



1X2) = ^ 
1X2/ AN 



(10) 



where (Si) is the GS expectation value of the i th - 
componcnt of the total spin of the system for the 
given choice of the quantization axis. On the basis 
of our discussion of the giant-spin model, (Si) can 



Jij 



(AN) 2 E^ 

ij 



(M% + M, 
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ij 



M 12 



(13) 



The tensor 7^ can now be diagonalized by a unitary 
transformation and A2 becomes 

A 2 = a + A((S' x )f + B((S' y )f + C((S' Z )) 2 (14) 
= a + A((S' X ) 2 ) + B((S' y ) 2 ) + C((S' z f) (15) 

where A, B,C are the eigenvalues of 7^ and the <S| 
are the three spin components rotated along its three 
principal axis (Eq. (fT5j) follows from Eq. (fT"4"|) thanks 
to the properties of spin coherent states). 

This expression for A 2 is exactly the expectation 
value (S, n| TL \S, n) of the quantum spin Hamilto- 



H = a + A(S' x y+B(S 



K) 2 



C(S. 



(16) 



in the spin coherent state \S, n). Eq. (JTHJ) is equiva- 
lent to Eq. © up to an irrelevant constant. 

The perturbative method described here works 
well for systems with a large HOMO-LUMO gap. 
However, for systems that have nearly degenerate 
and not fully occupied HOMO levels, which of- 
ten is the case for charged molecules, the pertur- 
bative approach breaks down, since some of the 
energy denominators in Eq. @ vanish. To avoid 
this problem the magnetic anisotropy can alter- 
natively be calculated by an exact diagonalization 
method. In this approach, the solutions of the one- 
particle Schrodinger equation in the SCF approxi- 
mation (which does not include SOI), are used to 
construct a finite matrix representation of the SOI, 
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which is then diagonalized exactly. The matrix is 
then diagonalized subject to the constraint that the 
resulting spin is aligned along a give choice of the 
quantization axis. The resulting single-particle so- 
lutions {e' k ,\ijj' k } = Y,<r WkJlXa)} are used to com- 
pute the trace of the system as a function of direction 
of the quantization axis (or average direction of the 
magnetization). In Ref3£, a discussion of the rela- 
tionship between the second-order variation in the 
trace and the self-consistent second-order variation 
of total energy is presented. 

Once one has obtained the trace as a function of 
axis of quantization, one can use relatively standard 
techniques to decompose the trace into a spheri- 
cal harmonic representation and then determine the 
magnetic principal axes. Alternatively by choos- 
ing magnetic principal axes that are equivalent to 
those predicted from the second-order expressions, 
it is always possible to directly compare exact- 
diagonalization results with the second-order results. 
Using exact diagonalization, One can further extract 
parts of the fourth- and higher-order anisotropy 
terms as well. However, since self-consistency and 
other terms also affect the magnetic anisotropy at 
fourth-order and beyond, the exact diagonalization 
results are primarily used to determined whether 
the second-order results are expected to be stable 
and a good approximation to experiment. In cases 
where near degeneracies occur at the Fermi level, 
the second-order and exact-diagonalization results 
can be very different especially if the states near the 
Fermi level are coupled by the spin-orbit interaction. 
For such much more careful analysis of re- 

sults is needed and it is reasonable to expect that 
some degree of self-consistency with non-collincar 
capabilities will be needed. 

For electronic-structure methods, such as NRL- 
MOL, where the wavefunctions arc expanded in 
terms of atom-centered localized basis functions 
the second-order perturbative method allows one to 
further analyze the anisotropy Hamiltonian on an 
atom-by-atom basis. By expanding the Kohn-Sham 
orbitals (\(f>kcr)) in Eq. © in terms of the atom- 
centered basis, the second-order expression (Eq. (0) 
can be decomposed into a sum over four centers^ 
The super-diagonal terms (all center indices the 
same) can then be used to determine an anisotropy 
Hamiltonian associated with each atom. In Reb 2 ^, 
this decomposition has been used to verify the per- 
pendicular hard-axis alignment model in the C04 
easy-plane magnetic molecule. In Refill, Baruah's 
method was used to demonstrate that essentially all 
of the magnetism in Mni2-acetate was due to the 
outer eight S~2 crown Mn ions, and that the sum 
of the single-ion anisotropics was very close to the 



total anisotropy. Further the degree to which non- 
additivity occurred was explained by a canting of 
the atom projected anisotropy axes relative to the 
global anisotropy axis. 

To complete the discussion about second-order 
anisotropy Hamiltonians, derived either perturba- 
tively or via exact diagonalization, it is important 
to note a contribution to van Wullcn and cowork- 
ersi 28 ' 29 . Van Wullen noted that once a method is 
used to determine the spin-orbit energy as a function 
of axis of quantization that an additional quantum 
correction is needed to determine the parameter's D 
and E in the anisotropy Hamiltonian. For example, 
the M s = eigenstate is not aligned with an axis 
of quantization along the x-axis, the y-axis or any 
other axis. Therefore more care must be taken to 
determine D once the classical energy as a function 
of expectation value of S is known. Accounting for 
this correction changes the definition of D, as origi- 
nally derived by Pederson and Khanna, by a factor 
of (S +1/2)/ S. While this correction is small in the 
large S limit, it can be important for systems with 
lower spin. 



C. Computational details 

In this paper we use a self-consistent field approx- 
imation based on density- functional theory (DFT). 
A review of this approach in the study of molecular 
magnets can be found in Refi 2 ^. Here we remind 
a few key features that arc relevant for the present 
work. 

In a DFT calculation of a magnetic molecule, 
we obtain the total energy of the system for spe- 
cially prepared spin symmetry-breaking metastablc 
states. In many cases these are ferrimagnctic spin 
configurations, suggested by experiment. The en- 
ergies of these different metastable states can then 
be compared and the lowest-energy spin configura- 
tion determined. Alternatively, it is also possible to 
impose a fixed spin configuration, which in princi- 
ple would not remain stable after convergence. In 
all these symmetry-breaking calculations the state 
with a given spin configuration is represented by a 
single Slater determinant of occupied single-particle 
states, constructed in terms of self-consistent Kohn- 
Sham eigenvectors. In the absence of SOI, the Kohn- 
Sham wavefunctions have a well-defined spin charac- 
ter, majority or minority spin. Therefore, the single 
Slater determinant, representing a given spin con- 
figuration, is an cigenstate of the component of the 
total spin S in the direction of the quantization axis 
n, which is the magnetization direction. In general 
this single Slater determinant is not an eigenstate 
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of S 2 , but in many cases it will have a large over- 
lap with the eigenstate of S 2 with S equal to the 
eigenvalue of S ■ n. The GS total spin S of the 
molecule (in the absence of SOI) is taken to be equal 
to one half of the excess of majority spin electrons, 
S = AN/2 = (N maj - N min )/2, for the metastable 
spin configuration with the lowest total energy. The 
spin magnetic moment of the system in units of the 
Bohr magneton /ib is then /is = ATV/ib = 2SfiB- 
In the DFT study of SMMs the possible presence 
of fractional occupancy of some of the KS wavefunc- 
tions close to the Fermi energy might result in nonin- 
teger values of -/V ma j — A^ m i n . Typically this happens 
when the HOMO-LUMO gap is very small or van- 
ishing. We will encounter examples of this difficulty 
in the study of the charged states of {Fe^}. The ex- 
istence of this general problem was first discussed by 
Janak et aS^ in reference to near degeneracies be- 
tween 3d and 4s electrons in neutral isolated atoms. 
In Rcf^i a set of equations are derived which, while 
cumbersome to solve, allow one to variationally de- 
termine the electronic occupations that satisfy the 
conditions proposed by Janak in Ref^. 

DFT can be used to extract the parameters defin- 
ing the spin Hamiltonian that is supposed to de- 
scribe the exchange interaction between the mag- 
netic ions of the molecule. First of all DFT can be 
used to ascertain whether or not there is a localized 
spin at each magnetic ion, by calculating the total 
spin polarization inside a sphere centered about a 
given atom. For typical SMMs, including the one 
considered in this paper, while the magnetization 
density is not localized entirely on the magnetic ions, 
the assumption of a well-defined quantum spin often 
turns out to be quite reasonable. Once this is estab- 
lished, the calculation of the total energy for a few 
spin configurations permits the computation of the 
exchange constants of Eq. (JTJ . We will see an exam- 
ple of this in the next section. 

The DFT calculations discussed herein are per- 
formed with the Gaussian-orbital-bascd NRLMOL 
program.] 32 ! 33 All calculations employ the Perdew- 
Burke-Ernzcrho^ (PBE) generalized-gradient ap- 
proximation for the density functional. A large basis 
set is employed in which the exponents for the single 
Gaussians have been fully optimized for DFT calcu- 
lations. The NRLMOL code employs a variational 
mesh for numerically precise integration and an an- 
alytic solution of Poisson's equation. 

All-electron calculations are performed for all el- 
ements of the {Fe^} SMM except for the Au atoms 
that are used to construct the leads attached to the 
molecule, for which we have used pseudo potentials. 
All the electronic and magnetic properties are cal- 
culated using an optimized geometry. 



III. ELECTRONIC AND MAGNETIC 
PROPERTIES OF ISOLATED {Fe 4 } SMM 

The chemical composition of the molecule used in 
this work is FeiC^H^O-^^ The four Fe atoms in 
{Fen} SMM form an equilateral triangle, as shown 
in Fig. [T] The molecule has idealized D% symme- 
try with the C2 axis passing through the central 
atom and one of the peripheral atoms. Using first- 
principles methods we have calculated, in detail, the 
electronic and the magnetic properties of the {F e 4 } 
SMM. 




FIG. 1. (Color online) Ball-stick top view of an isolated 
{Fen} SMM. Red, green, blue and purple balls corre- 
spond to iron, carbon, hydrogen and oxygen atoms, re- 
spectively 

In the ground state, the central Fe atom is coupled 
anti-fcrromagnctically with three peripheral atoms, 
whereas the three peripheral atoms arc coupled fcr- 
romagnetically with each other, as shown in Fig. O 
Each of the four Fe atoms has spin Sf & = 5/2, thus 
the total spin of the ground state is S = 5. The 
magnetic interactions between these atoms can be 
described by the Heisenberg spin Hamiltonian (Eq. 
©), 



H = J(si • s 2 + Si • s 3 + si • s 4 ) 

+ J'(s 2 • s 3 + s 3 • s 4 + s 4 • s 2 ) . (17) 

The two exchange parameters J and J' can be 
written in terms of the expectation values of the 
Hamiltonian of Eq. (|T7|) . for different spin config- 
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FIG. 2. Exchange interaction constants between four Fe 
atoms in {Fe 4 } SMM. 



urations. 
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J ~^zr(,E duuu E uuuu ) ) 
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Here, E duuu . E dduu and £ umra are the energies of 
the molecule where the spin orientations at their 
respective atomic positions (1,2,3,4) are labeled as 
d = down or u = up. Using NRLMOL we have 
calculated the energies of different spin configura- 
tion and upon substitution in Eq. (|17j) . we obtain 
J = 9.94 meV and J' = 0.64 meV. DFT calcula- 
tions overestimates J and J' since estimated values 
from susceptibility measurements^ are 2.62 meV 
and 0.14 meV, respectively. However, we note that 
the ratio of these two parameters agrees quite well 
for both theory and experiment. These values of 
the exchange constants ensure that the GS of the 
spin Hamiltonian Eq. (jTHJ) has indeed a total spin 
S = 5/2, well separated from excited states charac- 
terized by other values of S. 

Using first-principles methods we have calculated 
the electronic and magnetic properties of {Fe^} 
SMM for the neutral (Q = 0) and two charged 
states, namely, the anion (Q = —1) and the cation 
(Q = +1). (We will refer to the value Q as the 
extra charge added to the system.) A summary of 
the results is shown in Table HI These results can 
be understood with the help of the structure of the 
single-particle energy levels around the Fermi level 
in the absence of SOI, plotted in Fig. [3] for the three 
charge states Q = 0, ±1. 

The neutral molecule has a stable 5 = 5/2 GS, as 
anticipated. The HOMO-LUMO gap of the neutral 
isolated Fe^ molecule is is about 0.85 eV, where both 
HOMO and LUMO levels are minority (down) spin 



TABLE I. Properties of the isolated {Fe 4 } SMM for the 
three different charge states. *Note that the energy gap 
reported for the cation refers to the energy difference 
between the half-filled HOMO and empty LUMO. Sec 
Fig. Etc). 



Charge 
state 


Spin magnetic 
moment fis(^B) 


HOMO-LUMO 
energy gap (eV) 


Anisotropy 
barrier (K) 


Q=0 


10.0 


0.85 


16.05 


Q=+l 


9.3 


0.80* 


53.42 


Q=-l 


9.0 


0.06 


1.88 



-0.79 
-0.82 



-1.71 
-1.73 



-3.63 
-3.69 



(a) Q=-l 



4.65 
4.68 



-7.24 
-7.28 



.58 

'.92 



I Q=0 



(c) Q=+l 



FIG. 3. Energy levels (eV) of different charge states of 
isolated the {Fe^} SMM without spin-orbit coupling. Ej 
represents the Fermi level. The numbers on the right of 
the HOMO and HOMO-1 levels in (c) are the fractional 
occupancies of the corresponding level. 



states, see Fig. EJb). Apart from a small swapping 
of two levels below the HOMO, the level structure 
for the Q = — 1 charge state can be obtained from 
the the energy levels of the neutral system simply 
by filling the neutral LUMO with a spin-down elec- 
tron, see Fig. (Ha). As a result the total spin of the 
anion is S = 9/2. Note also that the HOMO-LUMO 
gap of the anion is now reduced to 0.06 eV com- 
pared to the neutral molecule. The electronic states 
changes significantly for the Q = +1 charge (cation) 
state: the two doubly degenerate spin-up HOMO-1 
and spin-down HOMO states of the neutral molecule 
swap place, see Fig. [3]Jc). Furthermore, since there 
is now one fewer electron, the new HOMO is now 
half-filled. This implies that the total spin of the 
cation is again reduced with respect to the neutral 
molecule by 1/2, that is, S = 9/2. The fractional oc- 
cupancy of the HOMO plays important role in the 
enhancement of the magnetic anisotropy, discussed 
below. Our DFT calculations yield values of the to- 
tal spin S in agreement with the level structure of 
Fig- El In particular for the charged states Q = ±1, 
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even when the initial spin configuration is set to be 
S = 11/2, the system converges eventually to the 
value 5 = 9/2. 

For the neutral molecule, the magnetic anisotropy 
landscape is characterized by an easy axis in the di- 
rection perpendicular to the plane containing the 
four Fe atoms (the z-axis). As shown in Table HI 
we find that the anisotropy barrier for this case is 
about 16 K, which is in agreement with previous 
calculations^ (In Ref2£, the authors have used two 
different {Fe^} complexes. The molecular symme- 
try for one of these complexes is C2, whereas the 
the other one has -D3 symmetry. Our calculations 
agree well with the one that has D3 symmetry). 
Note that a well-defined energy gap between occu- 
pied and unoccupied states (regardless of the spin), 
ensures that the perturbative and exact calculation 
of the anisotropy coincide. 

The magnetic anisotropy of the Q = +1 charge 
state has also uniaxial character in the z-direction, 
with a barrier of about 53 K, significantly larger than 
that obtained for the neutral molecule. On the other 
hand for the Q = — 1 charge state the anisotropy is 
reduced to about 1.9 K, with an easy axis in the 
XV-plane of the Fe atoms. The large change in the 
anisotropy for the two charged states has a very dif- 
ferent origin for the two cases and can be understood 
in the following way. 

For the Q = — 1 case, the small gap between 
the like-spin HOMO and LUMO states might at 
first suggest a breakdown of the perturbative treat- 
ment. In fact, these two states are coupled only 
minimally by SOI. The most important coupling oc- 
curs between the spin-down HOMO and the spin- 
up LUMO+1 states. The energy difference between 
these two states is « 0.25 eV. This value and the 
corresponding energy denominator in Eq. ^) is large 
enough for perturbation theory to work (as a com- 
parison with the exact approach clearly shows) and, 
at the same time, small enough for this individual 
transition to completely determine the main features 
of the anisotropy landscape. In particular, it turns 
out that this term in Eq. ^ favors an easy axis 
along a direction in the XF-plane of the Fe atoms. 
Since this magnetization direction is unfavorable for 
other terms in Eq. © (which prefer the z-dircction), 
there are positive and negative contributions in the 
total energy difference for the two magnetization di- 
rections (calculated with Eq. (J7J), which in the end 
lead to a reduced anisotropy barrier. 

The large enhancement of the magnetic anisotropy 
barrier found for the Q = +1 state can also be un- 
derstood with the help of single-particle energy dia- 
gram shown in Fig. [3jc). We observe that the half- 
filled doubly-degenerate spin-up HOMO level lies 



just above a (close-to-100%) filled doubly degenerate 
spin-down HOMO-1 level. The term involving tran- 
sitions between these two occupied and unoccupied 
levels totally dominate Eq.®. In fact, the small- 
ness of the corresponding energy denominator (a few 
meV) renders the perturbative approach inadequate. 
This is a classical example of a quasi-degeneracy at 
the Fermi level, where the exact treatment of SOI 
is necessary. As it is often the case, the inclusion of 
SOI lifts the quasi-degencracy for a particular direc- 
tion of the magnetization, leading to a substantial 
decrease of the total energy for that direction and, 
consequently, to a large anisotropy energy barrier. 



IV. {Fe 4 } SMM ATTACHED TO METALLIC 
LEADS 

In this section we investigate how the electronic 
and magnetic properties of the {Fe^} SMM change 
when the molecule is attached to metallic leads, as 
in transport experiments. The system that we have 
in mind is a single-electron transistor device, where 
metallic nano-leads, separated by a nanogap created 
by either break junction or electric migration, are 
bridged by a molecule functionalized with convenient 
chemical ligands. In our theoretical modeling, we 
are forced to find a convenient finite representation 
of otherwise semi-infinite leads in the form of finite 
clusters. For the calculations reported in this paper 
we have chosen to model a metallic nano-lead with a 
finite cluster of 20 gold ( 79 Au) atoms, arranged in a 
special tetrahedral structure, which can be viewed as 
a fragment of the face-centered cubic lattice of bulk 
Au (see Fig. @|. This metal cluster, Au 2 o, has been 
previously investigated by by Li et al. in Ref. HH, 
where it was shown that 20 Au atoms arranged in 
this geometric configuration form a very stable sys- 
tem. Its very large HOMO-LUMO gap (1.77 eV) 
makes Au 2 o chemically very inert. What is also 
important is that its unique tetrahedral structure 
makes this cluster an ideal model for Au surfaces at 
the nanoscale. Our rational for using this cluster to 
model metal leads is that during the fabrication of 
the nanometer-spaced electrodes, via self-breaking 
by electro-migration for example, the ensuing Au 
nano-leads will relax into the most stable configu- 
ration which might be well described by tetrahedral 
Au2o- In a way, the tetrahedral AU20 is the best 
representation of bulk Au at the nanoscale. 

Secondly, since ultimately we would like to inves- 
tigate transport properties of this system in coulomb 
blockade regime, it is essential that coupling between 
the leads and molecule is weak. 

After constructing the two leads in the form of 
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two AU20 clusters, we have connected the molecule 
via phenyl groups, as shown in Fig. [4J The function- 
alization of the ligands of {Fe^} SMM by means of 
phenyl groups is well-known and suitable to attach 
the molecule to Au surface o 12 i 15 . 

A similar way of connecting the {Fe^} SMM to 
Au leads is employed in SET experiments in the 
Coulomb blockade regime^—, and it ensures that 
the electronic coupling between the molecule and 
electrodes is weak. To maintain proper bonding and 
ionic neutrality of the entire cluster we have further 
removed one hydrogen (H) atom from the molecule 
close to the contact point and have added it to the 
chain part of the cluster. We have considered two 
different ways of doing this. In the first case (here- 
after called Type-1 lead) we have added the H atom 
to the sulfur (S) atom near the Au cluster. In the 
second case (hereafter called Type-2 lead) we have 
added H to the carbon (C) atom near the contact 
point as shown in Fig. @] After connecting the leads 
with the molecule we have relaxed the entire sys- 
tem again. Typically we find that the system with 
Type-2 leads is more stable, that is, its energy is 
approximately 0.6 eV lower than the energy of the 
system with Type-1 leads. We will nevertheless re- 
port results for both cases unless otherwise specified. 





Type 2 Lead Added H »*V^* Energy is lower than Type 1 



FIG. 4. {Fa} SMM connected to Au 2 o leads. Two types 
of leads are used in these calculations. In Type-1 lead a 
H atom is added to the S atom near the gold lead (top 
figure). In Type-2 lead a H atom is added to a C atom 
in the phenyl group near the {Fa} molecule (bottom 
figure). 

A summary of the magnetic properties of the neu- 
tral molecule attached to A20 leads is shown in Ta- 
ble |TTJ We can see that the combined molecule plus 
leads system maintains a sizable HOMO-LUMO gap 
of about 0.75 eV. The first important result is that 
coupling the leads does not cause a change of the 
spin of the molecule, which remains equal to the 



value of the isolated neutral {Fe^}, 5 = 5. Sec- 
ondly, the leads have a very small effect also on the 
magnetic anisotropy of the system: the magnetic 
anisotropy landscape has still an easy axis along the 
same z-direction (see white arrow in Fig. [4j , with an 
energy barrier quite close to the 16 K of the isolated 
molecule. 



TABLE II. Properties of neutral {Fa} SMM attached 
to A20 leads, compared to the properties of the isolated 
molecule (first row). Type-1 and Type 2 (called in the 
paper also Type-1 lead and Type-2 lead) refer to the 
two different choices to place an Hydrogen atom to the 
ligand, See Fig. [4] 





Spin magnetic 


HOMO-LUMO 


Anisotropy 




moment /xg(^s) 


gap (eV) 


barrier (K) 


{Fe 4 } 


10 


0.85 


16.05 


{Fe 4 }+Type 1 


10 


0.87 


15.99 


{Fe 4 }+Type 2 


10 


0.57 


15.47 



-4.06 



-4.43 I ]] — (Phenyl 

-4.44 — ^ — cy 

-4.72 - f 

l(Au„) 
-5.01 £ — 

-5.02 



-3.91 f 

-4.03 - f - 

-4.14 J-J 



T 



it 



■4.75 

4.78 



-5.80 -tf-U- 

(a) 



(b) 



Au lead Au+Fe 4 +phenyl 



(c) 

Fe 4 +phenyl 



FIG. 5. Shifts in the energy levels (eV) of neutral {Fa} 
SMM when connected to two AU20 leads (of Type-2). 
(a) Isolated Au lead (b) Leads+{7 7 'e4}+phenyl and (c) 
{i ? e4}+phenyl. The labels Aud and Phenyl C p in (b) 
indicate that the main contribution to those levels comes 
from d levels of the Au leads and p levels of C in the 
phenyl ligands respectively. 

We can gain some insight about the robustness 
of the magnetic structure of {Fe^} SMM under the 
influence of metallic contacts by investigating the 
changes in the single-particle energy levels and level 
alignment and occurring when the leads are con- 
nected to the molecule. We have calculated sep- 
arately the energy levels of the isolated Au lead(s) 
and the {F e 4 } + phenyl group, along with the levels 
of the combined system [Au leads + {Fa} + phenyl 
group]. The results are shown in Fig. [5] The states 
at and near the Fermi level of the two subsystems are 
dominated by the d-levcls of Au atoms and the p lev- 
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(a)Homo-ncutral 




(b)Lumo-ncutral 
FIG. 6. HOMO and LUMO of the neutral molecule. 



els of the C atoms of the phenyl group. Thus, when 
the two systems are combined, the charge transfer 
taking place to align the Fermi energies of the two 
subsystems is restricted only within the contact re- 
gion, leaving the magnetic properties of the {F e^} 
inner core unchanged. Fig. [5jb) shows that the en- 
ergy levels of the combined system around at the 
Fermi level correspond primarily to states remote 
from the magnetic core. This is also supported by 
Fig. |51 where we plot the probability density for the 
HOMO and LUMO states of the {Fe A } + leads sys- 
tem. Both states have negligible contributions on Fe 
atoms or atoms immediately nearby to these. As we 
will see below, this implies that the magnetic prop- 
erties will remain unchanged even when extra charge 
is added to (or subtracted form) the system. 

Fig. E^b) shows that both the non-degenerate 
spin-down HOMO and the non-degenerate spin- 
down LUMO of the neutral system lie quite close 
in energy to degenerate levels (occupied and non- 
occupied respectively). Therefore we can expect 
that subtle energy-level swaps may occur when one 
electron is added to or subtracted from the system. 
As it is evident from Fig. [7J this is exactly what 
the calculations show. For the case of leads of Type 
2, the HOMO of the anion (Q = -1) is now a half- 
occupied doubly-degenerate spin-up level, lying very 
close to a non-degenerate spin-down LUMO. This 
leads to a GS spin S = 11/2, and to a spin mag- 
netic moment close to fj, s = 11/ie (Sec Tabic ILTI)) . 39 
A similar situation occurs for the Q = +1 charge 
state, which has a GS spin S = 11/2. We find 
that this state is however almost degenerate with 
another state with S = 9/2. For leads of Type- 
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-3.63— f- 
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-4.37- 
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-5.92 | ^. -0.15 
-6.47— — ^ »-E f 

-6.48 f- 

-6.50- 
-6.53- 



a) anion b) neutral 



c) cation 



FIG. 7. Energy levels (eV) for the three charge states 
Q = 0, ±1 for the {Fe 4 } SMM connected to two Au 20 
leads (of Type-2). (a) Anion, (b) Neutral (same as in 
Fig. [Sfb)). (c) Cation. The labels Au d and Phenyl C p 
in (b) indicate that the main contribution to those levels 
comes from d levels of the Au leads and p levels of C in 
the phenyl ligands respectively. 



1 (which are less stable), the spin configurations 
S = 11/2 and 5 = 9/2 are almost degenerate for 
both charged states, Q = ±1. Note that the spin 
magnetic moment of the Q~ -1 charge state is now 
closer to fi$ = 9^b- The quasi-degeneracy of two 
different spin configurations is a situation where the 
assumption of the existence of a well-defined giant- 
spin model may not be entirely adequate. 



TABLE III. Magnetic properties of the three charge 
states when the {Fa} SMM is attached to Au leads 
as in Fig. [4] 



Charge 


Spin magnetic 


Anisotropy 


state 


moment /is(^i_g) 


barrier (K) 




Type 1 


Type 2 


Type 1 


Type 2 


Q=0 


10.0 


10.0 


15.99 


15.47 


Q=+l 


9.0 


10.95 


17.73 


14.74 


Q=-l 


9.6 


10.65 


11.23 


16.97 



As shown in the Table iHll in contrast to the case 
of the isolated {F e A } SMM where the anisotropy 
of the charged states are significantly different from 
that of the neutral molecule, when the leads are at- 
tached to the molecule the anisotropy barrier of the 
charged states remains close to the value of the neu- 
tral system. We also note that magnetic properties 
of the charge states have some dependence on type 
of the lead attached to the molecule. 

As anticipated above, an explanation of this be- 
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FIG. 8. (Color online) Evaluation of the fraction of extra 
charge for the anion state (Q = —1, one extra electron 
added) with respect to the neutral state, contained in 
different regions of the {Fe^} molecule plus leads and 
phenyl groups. In (a) each yellow sphere surrounding 
the lead and the phenyl group contains about 40% of 
an electron charge. In (b) each yellow sphere, surround- 
ing only the Au lead, contains 21% of electron charge. 
Therefore, the amount of charge transferred to the leads 
is about 42% and to the phenyl groups is 38%. The rest 
of the extra charge ~ 20% is in the {Fe^} region. 



havior is already suggested by the energy diagram of 
Fig. [SJb) and the plots of Fig. [5] demonstrating that 
both the HOMO and the LUMO of the neutral sys- 
tem are states predominately localized around the 
Au leads and within the phenyl group respectively. 
Therefore, we expect that when we add or remove 
an electron from the system it largely resides in the 
lead and phenyl group, leaving the magnetic states 
in {Fen} molecule relatively unchanged. The easy 
axis, in all cases, points perpendicularly to the Fe4 
plane, as shown in Fig. |4l except for Q = — 1 charge 
state of Type-1 lead, which is in the plane. 

Further support to this picture is provided by cal- 
culating the real-space location of the extra charge 
when an electron is added or subtracted to the sys- 
tem. As an example, we consider here the case of 
the anion, where one electron is added to the system. 
Since part of this extra charge might end up in inter- 
stitial regions between atoms (this is the case for the 
extra charge on the Au leads), particular care must 
be taken in drawing conclusions based only on the 
atomic-position plot of the HOMO states, shown in 
Fig. \§\ which might miss this contribution. To cap- 



(c)LUMO 

FIG. 9. The two quasi-degenerate HOMOs and the 
LUMO of the anion charge state. Approximately 20% 
of the HOMO wave functions reside on the Au leads, 
primarily on the interstitial space between Au atoms. 
See Fig. [S] This contribution is not visible on the scale 
of this plot. 



ture the interstitial contribution, we draw instead 
a large sphere enclosing a given region of the sys- 
tem. NRLMOL is able to calculate the extra charge 
contained globally in that region, including intersti- 
tial contribution. By repeating the same calculation 
with different spheres centered at different locations, 
we can eventually determine the amount of extra 
charge in different relevant parts of the system. 

we consider a sphere (yellow color) 



In Fig. 



containing both the lead and the phenyl group 
linker. We find that the amount of extra charge con- 
tained in this region is 40% of one electronic charge. 
The remaining 20% is located on atoms in the near- 
est surrounding of the {Fe±} core. In Fig. [8(b)] the 
sphere encloses only the Au lead but no linker. For 
this case we find that each Au leads contains 21% 
of extra electronic charge. We conclude that when 
one electron is added to the system, a total of 42%. 
of the extra charge resides on the leads, 38% on the 
ligands and only 20% is around the magnetic core 
of {Fe^}. This 20% of added charge in not directly 
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on the Fe atoms and therefore does not change the 
magnetic properties of {Fe^} significantly. 
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-3.883 -3.849 -3.815 -3.781 -3.747 
Spin magnetic moment 



(a)Central Fe atom 




3.8 3.85 3.9 3.95 

Spin magnetic moment 



(b)Vertex Fe atoms 

FIG. 10. Local magnetic moments of Fe atoms. The neg- 
ative and the positive moments in the two figures imply 
that moments of central and vertex atoms are opposite 
to each other. 

Further evidence of this important conclusion is 
provided by the comparison of the calculated lo- 
cal spin magnetic moments of the Fe atoms for 
the isolated {Fe^} molecule and for molecule-plus- 
lead system, for different charge states. The re- 
sults are shown in Fig [TUJ We note from the fig- 
ure that for the isolated molecule the magnetic mo- 
ments change considerably for the charged states, 
whereas for molecule plus lead system, the corre- 
sponding change is very small. Clearly, when the 
molecule is attached to the Au leads, adding or re- 
moving one unit charge affects the magnetic states 
of the {Fe^} minimally, which is also why we do not 
see a large change in magnetic anisotropy for differ- 
ent charge states. 

We conclude this section with a few comments on 
the important issues of the nature of the charged 
states and the character of the coupling molecule- 
leads, as evinced by the DFT calculations. Firstly, 
we have seen that when one electron is added or 
subtracted to the system {Fe^} + (finite) leads, the 
extra charge is predominantly localized on the lig- 
ands (rj 40%) and on the leads ( w 40%). If we 
could increase the size of the leads, a larger fraction 
of the extra charge would be likely to spread on the 



metallic leads. Therefore one could argue that the 
charged states Q = ±1 investigated above arc not a 
fully adequate description of the charged states in- 
volved in the sequential tunneling processes taking 
place in a SET, where the additional charge should 
be essentially localized on the central island. Sec- 
ondly, and partly connected to this issue, the non- 
zero amplitude of the LUMO wavefunction of the 
neutral system (which is quite close to the HOMO 
wavefunction of the anion) seems to indicate that 
the ligands considered here do not behave as tun- 
nel barriers of a SET, but rather model an exam- 
ple of strong coupling between molecule and leads^ 
Both these features could be due to limitations of the 
DFT approach considered here, which tends to over- 
delocalize any added charge. Such drawbacks can 
possibly be improved by more refined DFT treat- 
ments, involving, for example, self-inter corrections. 
While we believe that these refinements are impor- 
tant and should be further investigated, the goal 
and strategy of the present paper is to simulate with 
DFT a realistic example of SMM attached to leads, 
being aware of these limitations. 



V. EFFECT OF AN EXTERNAL ELECTRIC 
POTENTIAL 

In SET devices the charge of the central island 
weakly coupled to metallic electrodes can be varied 
experimentally one by one by applying an external 
electric field via a third gate electrode, which over- 
comes the charging energy e 2 /C of the island. Here 
we investigate the effect of an external gate elec- 
trode, whose electric potential tends to confine the 
extra charge closer to the molecule. Note that in 
phenomcnological studies of SETs based on model 
Hamiltonians, a gate voltage only shifts the energies 
of the isolated "quantum dot" without affecting its 
wavef unctions. As we show below, in our case the 
gate voltage can be used to localize a wavefunction 
closer to the molecule and modify its coupling to the 
leads. The resulting charged states should be a bet- 
ter representation of the states involved in tunneling 
transport in SET when Coulomb blockade is lifted. 

We model the external potential by a simple Gaus- 
sian confining potential of the form 

JJ = y ae -a x (x-x ) 2 -a y (y-y ) 2 ~a z (z-z ) 2 ^ j-^g-j 

Here Vq is magnitude of the potential centered at 
(xo,2/Oj^o): which in our case is the position of the 
central Fe atom of the {Fe^}. The constants a's are 
the width of the potential along the corresponding 
directions and are chosen so that the potential drops 
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quickly at distances larger than {Fe^}. The sign 
of Vq determines whether the extra electron will be 
confined into or repelled from the {Fe^} molecule. 
Thus for the anion case a negative Vo will attract 
the electron whereas for the cation case a positive 
Vo will attract the "hole" inside the molecule. 

We start by looking at the effect of the gate volt- 
age on the anisotropy of the isolated {Fe^} SMM. 
We have first considered a gate voltage that depends 
only on the variable z. The resulting electric field 
points along the the z-axis, which is the easy axis of 
the molecule. 

From Fig [IT] we note that in the anion case a con- 
fining potential for the extra electron (Vo < 0) re- 
duces the anisotropy barrier; whereas repelling the 
extra charge away from the molecule increases the 
anisotropy. The neutral molecule displays an op- 
posite behavior as a function of Vq. In both cases 
the the behavior of the anisotropy is close to a linear 
function of Vo < 0. As expected, the variation of the 
barrier for the neutral molecule is limited, less than 
10 % for the largest applied voltage. The cation is 
special. We have seen that at zero voltage, the sys- 
tem has a large anisotropy barrier (see Table HJ. due 
to a quasi-degeneracy at the Fermi level. The exter- 
nal potential lifts this degeneracy and the anisotropy 
barrier decreases sharply for both signs of Vq. 
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FIG. 11. Variation of the magnetic anisotropy barrier as 
a function of a confining potential Vo applied along C3 
axis (perpendicular to the plane of the {Fa} triangle) 
of an isolated {Fe 4 } SMM. 

The confining potential that we have applied 
above does not break the C 3 symmetry of the sys- 
tem and hence the parameter E in Eq. ©, charac- 
terizing the transverse component of the magnetic 
anisotropy, is zero. However, the symmetry can be 
broken by applying an electric field along directions 
other than the easy axis. Table IIVI shows the effect 
of this broken symmetry on the anisotropy of the 
isolated {Fe 4 } SMM. 



TABLE IV. The effects of confining potential on neutral 
{Fa} SMM. A is the anisotropy barrier and D and E 
are the parameters of the Hamiltonian Eq. ([2]), all in 
units of K. 



V 
(eV) 


ot x 


a v 


a z 


ANISOTROPY 


A (K) 


D (K) 


E (K) 


0.0 











16.06 


-0.63 


0.00 


-5.0 


0.01 





0.01 


17.00 


-0.64 


-0.03 


-5.0 


0.01 








15.77 


-0.60 


-0.02 



It is evident from Tablc llVl that E is no longer zero 
if the electric field is applied along directions other 
than the easy axis. A non-zero E allows different 
eigenstates of z-componcnt of the giant spin to mix 
with each other, a transverse component in principle 
can cause quantum tunneling of the molecule giant 
spin. Thus, this method can be used as electric con- 
trol of magnetic properties. It can play a significant 
role in transport, for example by modifying spin se- 
lections rule and by opening alternative channels via 
quantum tunneling of the magnetization. 

We now discuss the effect of the applied gate volt- 
age when the {Fe^} molecule is attached to Au 
leads. In this case we have applied the field only 
along the easy axis of the molecule attached to the 
leads of Type 2, as shown in Fig. SJ We have seen in 
the previous section that since HOMO and LUMO 
states and states close in energy to these are primar- 
ily localized on the Au leads and phenyl linker, an 
added or removed electron leaves {Fe^} largely un- 
affected. But the presence of a confining potential 
(Vo < for electrons), applied only on {Fe^} part of 
the system, brings the states localized within {Fe^} 
SMM closer to LUMO levels. Thus, when an elec- 
tron is added to the system, the fraction of this extra 
charge that goes inside the molecule increases as we 
increase the confining potential. Similarly, when an 
electron is subtracted from the system, an applied 
positive gate voltage (Vo > 0), tends to localize a 
fraction of the positive extra charge (a hole) closer 
to the molecule. 

As an example, we consider the effect of a con- 
fining potential for the anion case (Q = —1, one 
electron added to the system). Fig. [T2l shows the 
change in fractional charge that enters into the 
{Fei} molecule as the strength of confining poten- 
tial is increased, and the corresponding change in 
magnetic anisotropy barrier of the system. Clearly, 
as the voltage is increased, more of the added elec- 
tron is pushed inside the molecule. As the charge 
fraction approaches unity, the anisotropy decreases 
and converges to the value obtained for the anionic 
state of the isolated {Fe 4 } SMM. 

It turns out that not only the anisotropy barrier 
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to modify the magnetic properties of the SMM. This 
in turn can have important effects on the tunneling 
transport of the device. 
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FIG. 12. (Color online) Fraction of the added charge 
(blue circles) confined inside {Fe^}, as a function the 
confining potential strength, and corresponding change 
in the magnetic anisotropy (black diamonds) of the leads 
plus molecule system in the anion state, Q = — 1. The 
external potential corresponds to an electric field along 
the z-direction. See Fig. [4] 



but other magnetic properties converge to the prop- 
erties of the anion state of the isolated {Fe^} as the 
extra charge, under the effect of the external po- 
tential, moves closer to the center of the molecule. 
The easy axis of the system, which in zero potential 
points perpendicular to the plane of the Fe4 triangle 
and towards the leads, eventually rotates into the 
plane of the Fe4 triangle, exactly as in the case of 
the anion state of the isolated {F e^} SMM. Sim- 
ilarly, as the added charge moves inside the inner 
magnetic core of {Fe^}, the total spin of the system 
is reduced from 5 = 11/2 to the value of the anion 
state of the isolated molecule, 5 = 9/2. 

Similar results are obtained for the cation. As we 
apply an increasingly positive voltage, a larger frac- 
tion of a (negative) electron charge is pushed outside 
{Fe^}, or equivalently, a larger fraction of (positive) 
hole is attracted inside the {Fe^}. As a result, the 
anisotropy barrier increases and it reaches a value of 
22.8 K for Vq = 5 eV, with more than half of the 
extra (positive) charge now inside {Fe±\. Similarly, 
the spin also switches from 5 = 11/2 at Vq = to 
5 = 9/2 at Vq = 5 eV. Again, this is consistent with 
both the spin and the anisotropy converging towards 
the corresponding values of the isolated cation state. 

A summary of the dependence of the anisotropy 
barrier as a function of the external potential for all 
three charge state in shown in Fig. Q21 While the 
anisotropy of neutral state displays a weak depen- 
dence on the field, the anisotropy of the two charged 
states is significantly affected. These calculations 
demonstrate that, for a SET with a {Fe^} SMM as 
a central island, by manipulating the position of the 
additional charge with a gate voltage, it is possible 
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FIG. 13. Magnetic anisotropy as a function of a con- 
fining field for the molecule plus leads (Type 2) system, 
for the three charge states. The field is applied along 
the easy axis, that is, perpendicular to the plane of the 
molecule. The anisotropy for the the three charge states 
for the isolated molecule is included as a comparison (see 
diamond, circle and triangle symbols). 



VI. COMPARISON WITH SET 
EXPERIMENTS 

Recent SET experiment s] 16 ' 18 have permitted the 
first measurements of the magnetic characteristics 
of the {Fe^} SMM weakly coupled to Au leads, 
for both neutral and charged states. Comparison 
with our theoretical results has be done with cau- 
tion, since important details (e.g., the type of linker 
used, see below) might differ in the two cases. In 
a first study Zyazin et aL— have studied quantum 
transport in the inelastic cotunneling regime. By 
measuring the zero-field splitting of magnetic exci- 
tations and their dependence on the magnetic field, 
it was possible, with help of the model Hamilto- 
nian of Eq. @, to extract values of the giant spin 
5 and the magnetic anisotropy barrier A = DS 2 
for three adjacent charge states, N (neutral), N + 1 
(anion) and N — 1 (cation). For the neutral state, 
the spin of the molecule was found to be equal to 
expected value 5jv = 5 and the anisotropy bar- 
rier to be consistent with the value in the bulk 
phase, A at = 1.4mcV = 16.24K. These results are 
in good agreement with our theoretical estimate for 
the Q = charge state, see Table UTTl 

For the charged states, the situation is more com- 
plicated. For the reduced molecule (anion state, 
one electron added) the experimental measurements 
gave 5jv+i = 11/2 for the spin and Ajv+i = 
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2.7meV = 31.30K for the anisotropy. For the oxi- 
dized molecule (cation state, one electron removed), 
the measurements gave Sn-i = 9/2 and An-i = 
1.8meV = 20.90 K respectively. Comparing these 
findings with the results of our calculations (see Ta- 
bic HIU)) we can see that, for a given choice of lead 
type (Type 1 or Type 2), the theoretical value of the 
spin is consistent with the experimental one for ei- 
ther cation or anion, but not for both. We can also 
conclude that the experiment typically finds larger 
values for the anisotropy barrier, for both reduced 
and oxidized states, than the values predicted by 
theory. 

Several reasons can be responsible for these dis- 
crepancies. First, the functionalization of the 
molecule used in the experiment is slightly different 
from the one used in the calculations. In fact, in the 
experiments two different types of ligands were used. 
In one case (labeled as sample A) the {Fe4} SMM 
was connected to the Au lead the phenyl group. In 
the second case, (labeled as sample B) the {Fe^} 
was connected via a thiol group, CgS. The coupling 
molecule-lead turns out to be weaker in sample A 
than in sample B4i As shown above, in all our calcu- 
lations we have used a different type of linker, which 
was combination of a phenyl group and a thiol group. 

Secondly, in experiment different charge states are 
achieved by adding or removing electrons to the cen- 
tral region of SET via a gate voltage. In the theo- 
retical calculations, the two relevant charged states 
are constructing by adding or removing an electron 
to a system consisting of the {Fe^} connected to fi- 
nite leads. The extra charge is allowed to relax in 
the self-consistent field, and it occupies regions away 
from the {Fe±}, which affects the magnetic proper- 
ties of the system. Indeed, confining on the SMM 
with an external gate modifies the anisotropy bar- 
rier. 

Third, the evaluation of the magnetic properties 
from experiment done in Ref. [TH relies on the use 
of the model Hamiltonian of Eq. ^ . The fitting of 
the experimental results maintains a degree of un- 
certainty and arbitrariness, and moreover it could 
be problematic in cases of level degeneracy at the 
Fermi level, not uncommon for charged states. In 
this case, we have seen the the giant-spin model of 
Eq. @ might become inadequate. 

Finally, the method of Ref. [l6| relies on the mea- 
surement of inelastic cotunneling excitations, which 
is quite sensitive to the coupling between molecule 
and leads, and therefore it is a procedure not im- 
mune of uncertainties. Indeed, in a more recent pa- 
per Burzuri et al*^ introduced a novel gate- voltage 
spectroscopy technique which permits the measure- 
ment of the anisotropy of an individual SMM in dif- 



ferent charge states by tracking the dependence of 
the charge degeneracy point as a function of mag- 
netic field. This method is much more sensitive 
and accurate than the method based on conven- 
tional transport spectroscopy employed in Ref. [ill 
The spin Hamiltonian provides a good fit of the 
data if Sn = 5, Ajy = 16. 2K for the neutral state, 
SV+i = 9/2 and Ajy + i = 16. OK for the reduced 
state and Sn-i = 11/2 and Aat_i = 16. 5K for the 
oxides state. Furthermore, the orientation of the 
easy axis is found to exhibit only small variations 
among different charge states. Although some de- 
tails might differ ^2, these results are quite consistent 
with the small variations in anisotropy magnitude 
and the unchanged orientation of the easy axis that 
we find for the three charge states in our theoretical 
analysis. See Table HTT1 

As we mentioned above, the small variation in 
the anisotropy for different charged states found in 
our calculations is related to the fact that any ex- 
tra charge added to the molecule + leads system 
tends not to reside directly on the magnetic atoms, 
but mainly on the ligands and on the Au leads. We 
pointed out that this could be, in part, an artificial 
effect due to the way we constructed charged states 
in our finite-size system and to the delocalizing char- 
acter of our DFT approach. On the other hand, it 
is interesting that our estimates for the magnetic 
anisotropy are essentially consistent with the exper- 
imental results of Ref. [jji, which are obtained ex- 
actly at charge degeneracy points. At these special 
values of the external gate voltage the energy of two 
adjacent charge states is the same. We can imag- 
ine that at the degeneracy point the extra charge 
can swap energy-free from the electrodes and the 
molecule and might be localized primarily in region 
where the molecule is connected to the leads. If 
this were the case, the charge distribution shown 
in Fig. [S] could be in fact a more realistic descrip- 
tion of the charge states than previously anticipated. 
We could also surmise that, exactly as it happens 
in our calculations when we apply an external con- 
fining potential, changing the gate voltage to move 
away from the degeneracy point and make a given 
charge state more stable could strongly affect the 
the magnetic anisotropy. Indeed the results of the 
the first experiment^, where the anisotropy was not 
extracted at degeneracy points but in the middle of 
a Coulomb blockade diamond, show a significantly 
enhanced anisotropy for charged states. The two ex- 
perimental results could simply indicate that in one 
case the extra charge is localized closer to the SMM 
than in the other, exactly as it happens in our the- 
oretical modeling. 
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VII. CONCLUSIONS 

In this paper we have studied the electronic and 
magnetic properties of a {Fe^} SMM in a single- 
electron transistor (SET) geometry, using DFT 
as implemented in NRLMOL. We have modeled 
the system by a {Fe^} functionalizcd with phenyl 
groups attached to two metals leads described by 
AU20 nanoclusters. Our calculations show that 
the magnetic structure of the neutral {Fe^} SMM, 
that is, its spin ordering and magnetic anisotropy, 
remains stable in the presence of metallic leads. 
Specifically the ground state spin is S — 5 and the 
anisotropy barrier is of the order of 16K, like for 
the isolated {Fen}. This result is ascribed to the 
fact that, when attaching leads to {Fe±}, any charge 
transfer between the molecule and the metal leads 
occurs primarily in the contact region and on the 
ligands, but does involve the magnetic core of the 
molecule. 

Based on the properties of the HOMO and LUMO 
of the neutral system, when an electron is added 
or subtracted to the molecule-lead system, we find 
that the added charge (Q = =pl) is primarily lo- 
cated on the ligands and on the leads. As a result, 
while the total spin of this finite system changes 
by AS* = ±1/2, the magnetic anisotropy displays 
small variations both in magnitude and orientation 
with respect to the neutral state. In contrast, the 
anisotropy of the anion and cation states of an iso- 
lated {Fe^} is quite different from the values of neu- 
tral molecule, since the added extra charge pene- 
trates the region of the Fe atoms. The theoretical 
study of charged states Q = ±1 for the molecule- 
leads system is technically challenging, due to oc- 
currence of small HOMO-LUMO gaps and conse- 
quent fractional occupancies of the states around 
the Fermi level. Furthermore, DFT tends to over- 
delocalize any added charge in the peripheral parts 
of the system. Nevertheless, the analysis of these 
states presented here sheds light on the properties 
of a {Fe^} SET when individual electrons are added 
or subtracted to the "quantum dot" by overcoming 
the Coulomb charging energy with a gate voltage. 

We have shown that an external electric potential, 
modeling a gate voltage, can be used to manipulate 
the charge on the molecule-leads system and with 
that the magnetic properties of the device. In par- 
ticular, for the two charged states Q = ±1 when 
the extra charge, under the effect of the potential, is 
progressively removed from the ligands-leads region 
into the magnetic core of the molecule, the magnetic 
properties converge to the properties of the anion 
and cation states of the isolated {Fe^}. This is an 
example of the electric control of magnetism of a 



SMM in a SET. The charged states of the molecule- 
leads system in the presence of external fields stud- 
ied in this paper can be used to construct the tran- 
sition matrix elements entering a quantum master 
equation describing tunneling transport in a SET. 
With the limitations inherent to the DFT approach 
mentioned above, these states incorporate charging 
effects for the SMM weakly coupled to metal leads. 

We have compared the results of our numer- 
ical calculations with the results of two recent 
experimental studies of tunneling transport in a 
{Fe^} three-terminal device in the Coulomb block- 
ade regim o 16 ! 18 . This comparison must be made 
with caution since some important details (e.g., the 
precise atomic and electronic structure of the lig- 
ands) are different and might explain some of the 
discrepancies between theory and experiment that 
we find. Nevertheless, one of these experiments^ 
finds that the anisotropy for the two charged states 
Q = ±1 displays only small variations in magni- 
tude and orientation from the corresponding values 
of the neutral states, in agreement with our theo- 
retical findings. Interestingly enough, the experi- 
mental values are extracted by tracking the depen- 
dence of the charge degeneracy point between two 
adjacent charged states function of the mag- 
netic field. Our numerical calculations show that 
the nearly-indcpendence of the magnetic anisotropy 
on the charged states is related to the position of the 
added charge being far away from the magnetic core 
of the molecule. Thus the agreement between theory 
ad experiment might indicate that for a {Fe^} SET 
the charge of a added (subtracted) electron close to a 
charge-degeneracy point is primarily located on the 
ligands and in the contact region with the leads. If 
correct, this would be an example in which a mag- 
netic property of SMM-based SET can provide in- 
formation on the electronic properties of the charged 
states. 

For molecule-lead systems with finite gaps, we ex- 
pect our results to provide accurate predictions of 
experiment. However for those cases where HOMO- 
LUMO gaps are very small and the electronic states 
at the Fermi level are partially occupied, further 
understanding will require variationally accounting 
for the electronic occupations along the lines sug- 
gested in Refill Another point of view is that such 
fractionally occupied solutions are also strongly af- 
fected by self-interaction corrections and that ac- 
counting for such corrections will often significantly 
decrease the possibility of fractionally occupied so- 
lutions. Self-interaction corrections are also likely 
to provide a more complete understanding of the 
nature of charged states investigate in this paper. 
Addressing spin-dependent conductance in macro- 
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molecular to meso-scale devices will require efficient 
solutions to these problems and renewed efforts at 
extracting quantitative model Hamiltonians for such 
systems. 
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The fact that in Table Hill the spin magnetic moment 
of the Q = — 1 charge state is not exactly equal to an 
integer is due to the fractional occupancy (approxi- 
mately 15%) of the spin-down LUMO in our calcula- 
tions, carried out with a finite smearing of the Fermi- 
Dirac distribution. 

Note however that the amplitude of the LUMO wave- 
function is quite small in the central region. Therefore, 
an electron tunneling in from one of the leads would 



still find a bottleneck when tunneling out to the other 
lead. 

Note that in Ref . liH the results reported for the anion 
and cation were obtained for sample A and B respec- 
tively. 

In general, in Ref. [l8| it is found that upon reduction, 
either from N — > N + 1 or from N — 1 — s> N, the spin 
S always decreases and the anisotropy parameter D, 
defined via A = DS 2 , always increases. The results 
of our calculations (see Table [TTT|) show that both the 
anion and the cation have preferably 5* = 11/2 for a 
lead of Type-2, whereas S = 9/2 for a lead of Type- 
1. However, states with swapped spin configurations 
S —11/2 S = 9/2 are quite close in energy for 
both charged states. 



